Continuous- Variable Spatial Entanglement for Bright Optical Beams 
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A light beam is said to be position squeezed if its position can be determined to an accuracy 
beyond the standard quantum limit. We identify the position and momentum observables for bright 
optical beams and show that position and momentum entanglement can be generated by interfering 
two position, or momentum, squeezed beams on a beam splitter. The position and momentum 
measurements of these beams can be performed using a homodyne detector with local oscillator of 
an appropriate transverse beam profile. We compare this form of spatial entanglement with split 
detection-based spatial entanglement. 



I. INTRODUCTION 

The concept of entanglement was first proposed by 
Einstein, Podolsky and Rosen in a seminal paper in 1935 
PJ. The original Einstein-Podolsky-Rosen (EPR) entan- 
glement, as discussed in the paper, involved the position 
and momentum of a pair of particles. In this article, we 
draw an analogy between the original EPR entanglement 
and the position and momentum (x-p) entanglement of 
bright optical beams. 

Entanglement has been reported in various manifesta- 
tions. For continuous wave (CW) optical beams, these 
include, quadrature 0, Q an d polarisation Q entangle- 
ment. Spatial forms of entanglement, although well stud- 
ied in the single photon regime, have not been studied sig- 
nificantly in the continuous wave regime. Such forms of 
entanglement are interesting as they span a potentiall; 
infinite Hilbert space. Spatial EPR entanglement 
has wide-ranging applications from two-photon quan tum 
imaging |12lll3| to holographic teleportation |l llll5| and 
interferometric faint phase object quantum imaging |l6j| . 

Current studies are focused on x-p entanglement for 
the few photons regime. Howell et al. || observed 
near and far-field quantum correlation, corresponding to 
the position and momentum observables of photon pairs. 
Gatti et al. 6] have also discussed the spatial EPR as- 
pects in the photons pairs emitted from an optical para- 
metric oscillator below threshold. Other forms of spatial 
entanglement which are related to image correlation have 
also been investigated. A scheme to produce spatially en- 
tangled images between the signal and idler fields from an 
opticalparametric amplifier has been proposed by Gatti 
et al. Q ■ Their work was extended to the macroscopic 
domain by observing the spatial correlation between the 
detected signal and idler intensities, generated via the 
parametric down conversion process |lfj. 
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Our proposal considers the possibility of entangling 
the position and momentum of a free propagating beam 
of light, as opposed to the entanglement of local ar- 
eas of images, considered in previous proposals. Our 
scheme is based on the concept of position squeezed 
beams where we have shown that we have to squeeze the 
transverse mode corresponding to the first order deriva- 
tive of the mean field in order to generate the posi- 
tion squeezed beam |l7j . Similarly to the generation 
of quadrature entangled beams, the position squeezed 
beams are combined on a beam-splitter to generate x-p 
entangled beams. We also propose to generate spatial en- 
tanglement for split detection, utili sing sp atial squeezed 
beams reported by Treps et al. [ljl, |2(J, LLlJ ■ This form of 
spatial entanglement has applications in quantum imag- 
ing systems. 

We first define the position and momentum of an op- 
tical beam by performing a multi-modal decomposition 
on a displaced and tilted beam, respectively. We consider 
the case of a TEMoo beam and show that the correspond- 
ing position and momentum observables are conjugate 
observables which obey the Heisenberg commutation re- 
lation. We then propose a scheme to produce x-p entan- 
glement for TEMoo optical beams. Finally, we consider 
spatial squeezed beams for split detectors and show that 
it is also possible to generate spatial entanglement with 
such beams. 



II. POSITION-MOMENTUM ENTANGLEMENT 

A. Definitions - Classical Treatment 

Let us consider an optical beam with x- and y- sym- 
metric transverse intensity profile propagating along the 
z-axis. Since the axes of symmetries remain well de- 
fined during propagation, we can relate the beam po- 
sition relative to these axes. To simplify our analysis 
we henceforth assume without loss of generality, a one- 
dimensional beam displacement, d, from the reference 
x-axis (see Fig.^a)). We denote the electric field profile 
of the beam by E(x). For a displaced beam, the electric 
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FIG. 1: (a) Beam displacement d, and (b) beam tilt by angle 
6, from a reference axis. 



field profile is given by 



E d {x) = E(x) + d 



dE(x) d 2 d 2 E(x) 



dx 



dx 2 



(1) 



In the regime where displacement is much smaller than 
the beam size, we can utilise the linearised approxima- 
tion where only the zeroth and first order terms are sig- 
nificant. We see from this expression that the zeroth 
order term is not dependent on d, and that the displace- 
ment is directly proportional to the derivative of the field 
amplitude dE{x)/dx [13 • 

The transverse beam momentum p on the other hand, 
can be obtained from the transverse component of the 
wave-number of the beam, p — fcsin#, where k = 2ir/\ 
and the beam tilt is 9. This beam tilt is defined with 
respect to a pivot point at the beam waist, as shown in 
Fig.db). 

The electric field profile for a tilted beam with untilted 
electric field profile E(x) and wavelength A is given by 



Eg(x) = exp 



i2irx suit 
A 



Etxcost 



(2) 



We can again simplify Eq. @ by taking the zeroth and 
first order Taylor expansion terms to get a transverse 
beam momentum of p ~ k6. In the case of small dis- 
placement or tilt, we therefore obtain a pair of equations 



E d {x) 
E p (x) 



E(x) 
E(x) 



d 

dx 
p ixE(x) 



(3) 
(4) 



Eqs. © and give the field parameters that relate to 
the displacement and tilt of a beam. For freely propagat- 
ing optical modes, the Fourier transform of the derivative 
of the electric field, J-(dE(x)/dx), is of the form ixE(x). 
Hence, the Fourier transform of displacement is tilt. 

In the case of a single photon, the position and mo- 
mentum arc defined by considering the spatial proba- 
bility density of the photon, given by \E(x)\ 2 /I, where 
I = j \E(x)\ 2 dx is the normalisation factor. The mean 
position obtained from an ensemble of measurements on 
single photons is then given by (x) = y J x\E(x)\ 2 dx. 
The momentum of the photon is defined by the spatial 



probability density of the photon in the far- field, or equiv- 
alently by taking the Fourier transform of (x). These 
definitions are consistent with our definitions of position 
and momentum for bright optical modes. 



B. 



TEMp, Basis 



In theory, spatial entanglement can be generated for 
fields with any arbitrary transverse mode-shape. How- 
ever, as with other forms of continuous- variable entangle- 
ment, the efficacy of protocols to generate entanglement 
is highest if the initial states are minimum uncertainty. 
For position and momentum variables, the minimum un- 
certainty states are those which satisfy the Heisenberg 
uncertainty relation A 2 xA 2 p > %, in the equality. This 
equality is only satisfied by states with Gaussian trans- 
verse distributions [Lsl ] , therefore we limit our analysis to 
that of TEMoo modes. 

A field of frequency to can be represented by the posi- 
tive frequency part of the mean electric field £ + e luJt . We 
are interested in the transverse information of the beam 
described fully by the slowly varying field envelope £ + . 
We express this field in terms of the TEM pg modes. For 
a measurement performed in an exposure time T, the 
mean field for a displaced TEMoo beam can be written 



^)=i^^{u o{ x) + ^ { x)) (5) 

where the first term indicates that the power of the dis- 
placed beam is in the TEMoo mode while the second 
term gives the displacement signal contained in the am- 
plitude of the TEMio mode component. The correspond- 
ing mean field for a tilted TEMoo beam can be written 
as 



£ p +(x) 



TluJ 



iw p 



) + ^ Ul (x) (6) 



2e cT \ " ' 2 

where the second term describes the beam momentum 
signal, contained in the ir/2 phase-shifted TEMio mode 
component. 

C. Definitions - Quantum Treatment 

We now introduce a quantum mechanical representa- 
tion of the beam by taking into account the quantum 
noise of optical modes. We can write the positive fre- 
quency part of the electric field operator in terms of pho- 
ton annihilation operators a. The field operator is given 
by 



Tiuj 



2e cT ^ 

n— 



a n u n (x) 



(7) 



where u n (x) are the transverse beam amplitude functions 
for the TEM pg modes and d n are the corresponding an- 
nihilation operators. d n is normally written in the form 



3 



of a n = (a n ) + 5a n , where (a n ) describes the coherent 
amplitude part and Sa n is the quantum noise operator. 

In the small displacement and tilt regime, the TEMoo 
and TEMio modes are the most significant [I?) , with the 
TEMio mode contributing to the displacement and tilt 
signals. We can rewrite the electric field operator for 
mean number of photons N as 



St = 



sx+ + isx: 



2e cT \ 



E 

n=l 



iX; 



Ui(x) 



(8) 



where the annihilation operator is now written in terms of 
the amplitude and phase X~ quadrature operators. 

The displacement and tilt of a TEMoo beam is given 
by the amplitude and phase of the u\{x) mode in Eqs. (JSJ 
and ©, respectively. We can, therefore, write the beam 
position and momentum operators as 



wo 



2VN 
1 



x: 



WqVN 



x: 



(9) 
(10) 



D. Commutation Relation 

Two observables corresponding to the position and mo- 
mentum of a TEMqo beam have been defined. We have 
shown that the position and momentum observables cor- 
respond to near- and far-field measurements of the beam, 
respectively. Hence, we expect from Eqs. © and (110(1 . 
that the position and momentum observables do not com- 
mute. Indeed, the commutation relation between the two 
quadratures of the TEMio mode is = 2i. This 

leads to the commutation relation between the position 
and momentum observables of an optical beam with iV 
photons 



an equivalent dependence of v N. This validates the fac- 
tor of N in the Heisenberg inequality and the commuta- 
tion relation for a CW optical beam. 



E. Entanglement Scheme 

We have shown that the position and momentum ob- 
servables of CW TEMoo optical beams satisfy the Heisen- 
berg commutation relation. Consequently, EPR entan- 
glement for the position and momentum of TEMoo beams 
is possible. Experimentally, the usual quadrature entan- 
glement is generated by mixing two amplitude squeezed 
beams on a 50:50 beamsplitter. The same idea can be 
applied to generate EPR x-p entanglement, by using po- 
sition squeezed beams [17|. Our scheme to produce x-p 
entanglement between two CW TEMoo optical beams is 
shown in Figure El The position squeezed beams in (a) 




FIG. 2: Scheme for generating position-momentum entangle- 
ment for continuous wave TEMoo optical beams. OPA: opti- 
cal parametric amplifiers for the generation of squeezed light, 
BS: 50:50 beam-splitter, HD: homodyne detectors, LO: local 
oscillator beams, <f>: phase of local oscillator beam. 
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(11) 



This commutation relation is similar to the position- 
momentum commutation relation for a single photon, 
aside from the 1/N factor. The 1/N factor is related to 
the precision with which one can measure beam position 
and momentum. Rewriting the Heisenberg inequality us- 
ing the commutation relation, gives 



A 2 £A 2 P>lA7. 



(12) 



The position measurement of a coherent optical beam 
gives a signal which scales with TV. The associated quan- 
tum noise of the position measurement scales with >/~N. 
Hence the positioning sensitivity of a coherent beam 
scales as \/~N frUITi^ . The same consideration applied to 
the sensitivity of beam momentum measurement shows 



and (b) are generated via the in-phase combination of a 
vacuum squeezed TEMio beam with a coherent TEMoo 
beam. Such beam combination can be achieved experi- 
mentally, for example using an optical cavity or a beam- 
splitter [21|. The result of the combination is a position 
squeezed beam. To generate entanglement, we consider 
beams with zero mean position and momentum, but we 
are interested in the quantum noise of the position and 
momentum of the beam. With this assumption, the elec- 
tric field operators for the two input beams at the beam- 
splitter are given by 



c 2 



huj 

2e cT 



(y/Nu {x) + 



5a r , 



2e cT V 



Nuo(x) 



n=0 



Sb n u n (x) 



(13) 



(14) 
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where in both equations, the first bracketed term de- 
scribes the coherent amplitude of the TEM 00 beam. The 
second bracketed terms describe the quantum fluctua- 
tions present in all modes. For position squeezed states, 
only the TEMio mode is occupied by a vacuum squeezed 
mode. All other modes are occupied by vacuum fluctua- 
tions. It is also assumed that the number of photons in 
each of the two beams, during the measurement window, 
is equal to N. The two position squeezed beams (1,2) are 
combined in-phase on a 50:50 beam-splitter (BS) in (c). 

The usual input-output relations of a beam-splitter 
apply. The electric field operators describing the two 
output fields from the beam-splitter are given by £3 — 
(£1 + £-i)j\f2 and £4 — (£1 — £-i)j\f2. To demonstrate 
the existence of entanglement, we seek for quantum corre- 
lation and anti-correlation between the position and mo- 
mentum quantum noise operators. The position operator 
corresponding to beams 3 and 4 are given respectively by 



5x 3 = p— Ux++SX- 
8x4 = ^— (sx+ -5Xr 



(15) 



(16) 



The momentum operator corresponding to the photo- 
current difference for beams 3 and 4 are given by 



5p3 = 



Sp 4 = 



1 



Wo 



SX-+6X+) 
1 ( t - 2 



V2VN 



w 





8X--6X+) 



-^=(Sp a %Sx b 



(17) 



(18) 



In our case where the two input beams are position 
squeezed, the sign difference between the position noise 
operators in Eqs. H15fl and (|16fl as well as that between 
the momentum noise operators in Eqs. I|17l) and (|18fl are 
signatures of correlation and anti-correlation for 5x and 
5 p. 



F. Inseparability Criterion 

Many criterions exist to characterise entanglement, for 
example the inseparability criterion |22| and the EPR 
criterion [2j|. We have adopted the inseparability cri- 
terion to characterise position-momentum entanglement. 
For states with Gaussian noise statistics, Duan et al. |22| 
have shown that the inseparability criterion is a necessary 
and sufficient criterion for EPR entanglement. 



In the case where two beams are perfectly interchange- 
able and have symmetrical fluctuations in the amplitude 
and phase quadratures, the inseparability criterion has 
been gene ralised and normalised to a product form given 

by sin 111 mill Ei 

for any pair of conjugate observables x and p, and a pair 
of beams denoted by the subscripts 3 and 4. For states 
which are inseparable, T(5x,5p) < 1. By using observ- 
ables 5x and 5p from Eqs. (JTSJ, JTHJl, J32J and JTHJ) as 
well as the commutation relation of Eq. (fill) the insepa- 
rability criterion for beams 3 and 4 is given by 



I{5x, dp) 



16A 2 



+ ^2\ 



((sx+y}((5x+y) 



< 1 



(20) 



where we have assumed that the TEM10 modes of beams 
1 and 2 are amplitude squeezed (i.e. ((SX^) 2 ) < 1 and 
<(5X+) 2 )<1). 

Thus we have demonstrated that CV EPR entangle- 
ment between the position and momentum observables 
of two CW beams can be achieved. 



G. Detection Scheme 

Ref. has shown that the optimum small displace- 
ment measurement is homodyne detection with a TEM10 
local oscillator beam (see Fig. 121(d)). When the input 
beam is centred with respect to the TEM10 local oscil- 
lator beam, no power is contained in the TEM10 mode. 
Due to the orthogonality of Hermite-Gauss modes, the 
TEM10 local oscillator only detects the TEM10 vacuum 
noise component. As the input beam is displaced, power 
is coupled into the TEM 10 mode. This coupled power 
interferes with the TEM10 local oscillator beam, causing 
a change in photo-current observed at the output of the 
homodyne detector. Thus the difference photo-current of 
the TEM10 homodyne detector is given by [l7j 



2VNy/N ] 

w 



LO 



(21) 



where -/Vlo and N are the total number of photons in the 
local oscillator and displaced beams, respectively, with 
A^lo 3* N. The linearised approximation is utilised, 
where second order terms in 5a are neglected since N 3> 
I (Sa^ ) I for all n. 

In order to measure momentum, one could use a lens to 
Fourier transform to the far-field plane, where the beam 
is then measured using the TEM10 homodyning scheme. 
However, we have shown that the the position and mo- 
mentum of a TEMqo beam differs by the phase of the 
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TEMio mode component. Indeed for a tilted TEMoo 
beam, the TEMio mode component is 7r/2 phase shifted 
relative to the TEMoo mode component. Consequently 
the phase quadrature of the TEMio mode has to be in- 
terrogated. This can be achieved by utilising a TEMio 
local oscillator beam with a it/ 2 phase difference relative 
to the TEMio mode component of the TEM o beam. The 
resulting photo-current difference between the two homo- 
dyning detectors, for TVlo ^> N, is given by 

h p _ = woVn^/K^P (22) 

III. SPATIAL ENTANGLEMENT FOR SPLIT 
DETECTION 

The entanglement presented in the previous section 
is analogous to x-p entanglement in the single photon 
regime. However, the choice of the mean field mode is 
restricted to the TEMoo mode. This limits the richness of 
a spatial variable and thus excludes the possibility of gen- 
crating an infinite Hilbert space. To exploit the proper- 
ties of spatial variables, we now consider more traditional 
forms of spatial squeezing. Consequently, we study the 
possibility of generating spatial entanglement for array 
detection devices, based on spatial squeezed beams. 

A. Spatial Squeezing 

Spatial squeezing was first introduced by Kolobov [25| . 
The generation of spatial squeezed beams for split and ar- 
ray detectors was experimentally demonstrated by Treps 
et al. one-dimensional spatial squeezed 

beam has a spatially ordered distribution, where there 
exists correlation between the photon numbers in both 
transverse halves of the beam. A displacement signal ap- 
plied to this beam can thus be measured to beyond the 
QNL. 

We consider a beam of normalised transverse ampli- 
tude function vq(x) incident on a split detector. The 
noise of split detection has been shown to be due to the 
Sipped mode |24|, given by 

,,W_J«dW forxX) 
Vl[X) ~ \ -v (x) fovx<0 

When the field is centred at the split-detector, such 
that the mean value of the measurement is zero, the 
flipped mode is thus orthogonal to the mean field mode. 
In this instance, modes Vi(x) (for i > 1) can be derived 
to complete the modal basis. The electric field operator 
written in this new modal basis is given by 



V 2e cT V ^ ) 

where the first term describes the coherent excitation of 
the beam in the vq(x) mode and N is the total number 



of photons in the beam. It has been shown that the cor- 
responding photon number difference operator for split 
detection is given by 

n W = </N8X+ (24) 

The beam is spatially squeezed if the state of the 
flipped mode is vacuum squeezed and in phase with the 
mean field mode (see Fig. 01(a) and (b)). 



B. Spatial Homodyne 

Since split detection is commonly used as a detection 
device for beam position, one would naturally consider 
taking the Fourier transform of a spatial squeezed beam 
to obtain the conjugate observable for the beam. How- 
ever, we have shown that split detection does not cor- 
respond exactly to beam position measurement. Thus 
the Fourier plane of the spatial squeezed beam does not 
provide the conjugate observable. More practically, the 
flipped mode is not mode-shape invariant under Fourier 
transformation. In the far-field, each odd-ordered mode 
component of the flipped mode obtains a (2n + l)7r Gouy 
phase difference, compared to the near-field. Thus the 
mode-shape in the far-field is no longer a flipped mode. 
Consequently, far and near-field measurements of a spa- 
tial squeezed beam will not give the conjugate observ- 
ables. 

However, we can find the conjugate observables of a 
spatial squeezed beam by drawing an analogy to stan- 
dard homodyne detection. In split detection, the equiv- 
alent local oscillator mode is the mean field Vq(x) mode. 
The mode under interrogation by the split detector is the 
flipped mode v\ (x) . In the case of homodyne detection, 
the phase of the local oscillator beam is varied to mea- 
sure the conjugate observables (i.e. amplitude and phase 
quadratures) of the input beam. Adapting this concept 
to the split detector, the conjugate observables for the 
spatial squeezed beam is thus the amplitude and phase 
quadratures of the flipped mode, while the mode-shape 
of the flipped mode remains unaltered. This is further 
verified upon inspection of Eq. l(2^|) . 

Our scheme to perform a phase measurement of the 
flipped mode is shown in Fig-El (d) . In our scheme we as- 
sume that the mean field is a TEMoo mode. Note that in 
principle, this analysis could be performed for any mode- 
shape. The coherent TEM o mode component provides a 
phase reference for the flipped mode, analogous to that of 
a local oscillator beam in homodyne detection. Thus the 
phase quadrature of the flipped mode can be accessed by 
applying a it/2 phase shift between the the TEMoo mode 
and the flipped mode noise component. Experimentally, 
this is achievable using an optical cavity. When the cav- 
ity is non-resonant for the vq(x) and v\(x) modes it will 
reflect off the two modes, in phase, onto the split de- 
tector. This will give a measurement of the amplitude 
quadrature of the flipped mode. However, the cavity can 
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be tuned to be partially resonant on the vq (x) mode while 
reflecting the flipped mode. A tt/2 phase difference can 
then be introduced between the reflected vq(x) and Vi(x) 
modes, giving a measurement of the phase quadrature of 
the flipped mode. The corresponding photon number dif- 
ference operator is 

nL" 3 = VN8X~ (25) 

which is the orthogonal quadrature of the spatial 
squeezed beam. The photon number operators corre- 
sponding to the two measurements in Eqs. (|24[1 and i|25|) 
are conjugate observables and satisfy the commutation 

relation [n^,nL _) ] = 2iN. 

It is important to realise that the number of photons 
N in Eqs. I|24[l and (|25l) are only approximately equal. 
This is due to the fact that partial power in the TEMoo 
mode is transmitted by the cavity, when the cavity is 
partially resonant on the TEMoo mode. Although it is 
possible to implement a scheme that conserves the total 
number of photons at detection (e.g. losslessly separating 
the mean field and flipped modes and recombining them 
with a phase difference) , we would like to emphasise that 
our scheme is more simple and intuitive, as well as being 
valid when N is large. 

C. Entanglement Scheme 

In order to generate spatial entanglement for split de- 
tection, two spatial squeezed beams labelled 1 and 2 
are combined on a 50:50 beam-splitter, as shown in Fig- 
ure 09(c). 




FIG. 3: Scheme for generating spatial entanglement for 
TEMoo continuous wave light beams. OPA: optical paramet- 
ric amplifiers for generating squeezed beams, PP: phase plates 
for producing flipped modes, BS: 50:50 beam-splitter. 

The electric field operators for the two input spatial 
squeezed beams at the beam-splitter are described in a 
form identical to that of Eq. (|23[1 . The annihilation op- 
erators of the electric field operators for input beams I 



and 2 are labelled by c„ and d n , respectively. By fol- 
lowing a similar procedure as before, the photon number 
difference operator for output beams 3 and 4 from the 
beam-splitter are calculated. 

For the amplitude quadrature measurement, the addi- 
tion of the difference photo-current between beams 3 and 
4 yields 

4- 5 + nt ] = + SX+ ) (26) 

For the phase quadrature measurement, the subtraction 
of the difference photo-current between beams 3 and 4 
gives 

4- ] - nt ] = VN(8X+ - 6X+ ) (27) 

To verify spatial entanglement, the inseparability cri- 
terion is utilised. The substitution of Eqs. (|2*6l . lt2*7|) 
and the commutation relation between the photon num- 
ber difference operators into the generalised form of the 
inseparability criterion gives 

< 1 (28) 

where V Co = ({SX+) 2 } and V d() = ((SX+) 2 ) are the 
variances for the flipped mode component of the spatial 
squeezed beams 1 and 2. The inseparability criterion is 
satisfied for amplitude squeezed flipped modes V Co < 1 
and Vdo < 1. 

We have proposed a scheme to generate spatial entan- 
glement for split detection, using spatial squeezed beams. 
Spatial squeezing has been defined for any linear mea- 
surement performed with an array detector |26| . Sim- 
ilarly, spatial entanglement corresponding to any linear 
measurement, can be obtained. For an infinite span ar- 
ray detector with infinitessimally small pixels, it is thus 
possible to generate multi-mode spatial entanglement, in- 
creasing the Hilbert space to being infinite-dimensional. 

IV. CONCLUSION 

We have identified the position and momentum of a 
TEMoo optical beam. By showing that x and p are con- 
jugate observables that satisfy the Heisenberg commu- 
tation relation, a continuous variable x-p entanglement 
scheme is proposed. This proposed entanglement, as con- 
sidered by EPR 0, was characterised using a generalised 
form of the inseparability criterion. 

We further explored a form of spatial entanglement 
which has applications in quantum imaging. The detec- 
tion schemes for quantum imaging are typically array de- 
tectors. In this article, we considered the split detector. 
We utilised the one-dimensional spatial squeezing work 
of Treps et al. and proposed a spatial homodyning 
scheme for the spatial squeezed beam. By identifying 
the conjugate observables for the spatial squeezed beam 
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as the amplitude and phase quadratures of the flipped 
mode, we showed that split detection-based spatial en- 
tanglement can be obtained. 
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